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Assignment of resonances in dissociative recombination of HD+ ions: high-resolution 
measurements compared with accurate computations 
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The collision-energy resolved rate coefficient for dissociative recombination of HD"'" ions in the 
vibrational ground state is measured using the photocathode electron target at the heavy-ion stor- 
age ring TSR. Rydberg resonances associated with ro-vibrational excitation of the HD"*" core are 
scanned as a function of the electron collision energy with an instrumental broadening below 1 meV 
in the low-energy limit. The measurement is compared to calculations using multichannel quantum 
defect theory, accounting for rotational structure and interactions and considering the six lowest 
rotational energy levels as initial ionic states. Using thermal equilibrium level populations at 300 
K to approximate the experimental conditions, close correspondence between calculated and mea- 
sured structures is found up to the first vibrational excitation threshold of the cations near 0.24 
eV. Detailed assignments, including naturally broadened and overlapping Rydberg resonances, are 
performed for all structures up to 0.024 eV. Resonances from purely rotational excitation of the 
ion core are found to have similar strengths as those involving vibrational excitation. A dominant 
low-energy resonance is assigned to contributions from excited rotational states only. The results 
indicate strong modifications in the energy dependence of the dissociative recombination rate coef- 
ficient through the rotational excitation of the parent ions, and underline the need for studies with 
rotationally cold species to obtain results reflecting low-temperature ionized media. 

PACS numbers: 34.80.Gs, 34.80.Ht, 34.80.Lx 



I. INTRODUCTION 

The dissociative recombination (DR) process [1] of di- 
atomic ions, together with its competing reactions - ro- 
vibrational excitation and de-excitation as well as, at 
higher energy, dissociative excitation - plays a decisive 
role in astrophysical [2] and planetary [3] ionized me- 
dia, fusion plasma in the divertor region [3] and in most 
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cold plasmas of technological interest [5]. The hydro- 
gen cation and its isotopologues represent benchmark 
systems frequently taken up in experimental and theo- 
retical studies. With the technique of merged electron 
and ion beams [B] , collision studies can be realized under 
quasi-monoenergetic conditions for detailed comparisons 
of experiment and theory. This yields collision-energy re- 
solved rate coefficients while, for sensitive comparison to 
theory, good control is needed over the ro-vibrational ex- 
citation of the incident cations. Following first results on 
HJ in single-pass merged beams |7], the infrared active 
isotopologue HD+ moved into the focus of experiments 
[5HTT] using the ion-storage ring technique [H], where 
fast ion beams are stored over many seconds and vibra- 
tional de-excitation of HD+ is completed within much 
less than a second [13] ■ On H J , storage-ring experiments 
gave access to vibrational de-excitation in electron col- 
lisions [T4l|T5]. Moreover, using special ion sources, the 
control of the initial states could be much improved in re- 



cent experiments 



on Hj -|- e . They resulted in 



accurate state-to-state, energy-resolved rate coefficients 



whose comparison to advanced theoretical studies [131 [TO] 
is underway. 

Studies on the system HD+ + e~ have reached a high 
level of accuracy. Storage-ring experiments were com- 
pared to detailed predictions 20t i21j from multichannel 
quantum defect theory (MQDT) that addressed in par- 
ticular the dependence of the rate coefficient on the ro- 
tational excitation of the ion. The energy-resolved DR 
rate coefficient [TT] and rotational effects in the resonant 
collision processes |22] were experimentally investigated. 
Moreover, experimental and theoretical results were com- 
pared for the related process of rotational de-excitation 
in low-energy electron collisions ^23, . 

The study presented here on the DR reaction 



HD^ 




(1) 



is motivated by the large improvement of the experimen- 
tal energy resolution achieved by using a photocathode 
electron beam j24j for merged-beams experiments |25j . 
in which electron temperatures below 1 meV are reached. 
As a function of the collision energy, scanned with high 
resolution, the DR rate coefficient displays a rich struc- 
ture of resonances, caused by resonant capture of the in- 
cident electron into Rydberg states on a ro-vibrationally 
excited HD+ core before the dissociation of the recom- 
bined system into the channels given by Eq. (fTj). These 
measurements will be presented together with a detailed 
comparison to state-by-state theoretical calculations of 
the DR rate coefficient by the MQDT method with the 
goal to assign the dominant resonance structures. The 
assignment provides a highly sensitive test for the qual- 
ity of the molecular data - potential energy surfaces and 
mutual couplings - and for the theoretical approach to 
account for the underlying mechanisms and interactions. 
The setup allowing our current measurements is de- 



scribed in Sec. Hll Section III presents our experimental 



results. The main steps in the MQDT computation of the 
cross section are reviewed in Sec. |IV[ and the procedure 
of resonance assignment is described in Sec.|VJ A detailed 
comparison between measurements and computations is 
performed in Sec. |VI[ 



II. 



EXPERIMENTAL SETUP 



The experiment was performed in the heavy ion stor- 
age ring TSR at the Max Planck Institute for Nu- 
clear Physics in Heidelberg, Germany. IID+ ions from 
a Penning source were accelerated in a radiofrequency 
quadrupole accelerator and a linear rf accelerator to a fi- 
nal energy 5.3 MeV, injected into the storage ring within 
30;US, and stored for up to 10 s before their replacement 
by a new injection. Along their closed orbit with a cir- 
cumference of 55.4 m, the residual gas pressure is below 
1 X 10"^'^ mbar. As in previous experiments [2 6) at this fa- 
cility, the ions were merged on each turn in the ring with 



collinear, velocity-matched or slightly velocity-detuned 
electron beams of ~1.5m length in two straight sections 
of the TSR: the electron cooler p7H25] , and the electron 
target section (ETS) [5D|. 

The electron cooler was used to compress the phase- 
space density of the injected ion beam within the ffi'st 
1-2 s after the injection, reaching ^Imni ion beam di- 
ameter, much smaller than the electron beam diameters 
of ~27 mm and ~19 mm in the electron cooler and in the 
ETS, respectively. The electron cooler with beam veloc- 
ity Ve^c locked the ion beam velocity to v^ — v^.c] the lab- 
oratory electron energy of this electron beam was near 
970 eV. It was generated by a thermionic cathode and 
magnetically expanded so that the temperature in the 
co-moving frame of the beam can be calculated [3S1 [51] 
to ksT^ = 16 meV. This temperature refers to the trans- 
verse direction perpendicular to the ion beam, while the 
longitudinal electron velocity differences correspond to 
a temperature found to be ksTu ~ 0.06 meV in inde- 
pendent measurements [31j . The electron density of the 
electron cooler amounted to ne,c=l-8xl0^ cm~'^ and a 
longitudinal magnetic field of 0.04 T was present in its 
interaction region with the ion beam. 

The electron beam in the ETS was used for the DR 
rate coefficient measurements at variable beam velocity 
Ve.t, which sets the electron-ion collision energy through 
the detuning energy defined as 



Ed = {me/2) {Vi - Ve,t) 



(2) 



where m^ is the electron mass. The electron-ion collision 
velocities are described by an anisotropic thermal dis- 
tribution [32] with the transverse and longitudinal elec- 
tron temperatures T± and Tj|, respectively. The elec- 
tron beam of the ETS was produced with a photocath- 
ode at liquid-nitrogen temperature [24] and, similar to 
the electron cooler, magnetic expansion was applied be- 
tween the cathode and the interaction region [30] . At 
matched beam velocities {ve,t = Vi = Ve,c, correspond- 
ing to Ed = 0) the electron density of the ETS was 
Tie^t — 3.2 X 10^ cm^"^; a longitudinal magnetic field of 
0.057 T was present in its interaction region. Electron 
temperatures of ksT^ — (0.6 ± 0.2) meV and fc^Tn near 
25(UeV were determined for the photocathode electron 
beam in the present measurement, observing energy de- 
pendences in the DR rate coefficient as described below. 
The neutral fragments emerging from the electron-ion 
reactions in the ETS travel with the ion beam velocity 
and leave the storage ring in forward direction towards 
a dedicated beamline in which a surface barrier detec- 
tor (SBD) with 48x48 mm^ active area is situated at 
about 11.5 m from the ETS center. The SBD has unit 
efficiency and yields signals proportional to the kinetic 
energy of the fragments impinging within the sampling 
time of 1 /is which is long compared to the arrival time 
difference of fragments from one DR event of ~1 ns, while 
short compared to the typical interval between two in- 
dependent events at rates of 1 kHz. Signals associated 
with the full mass of HD+ are associated with a DR 



event, whereas those corresponding to smaller masses are 
treated as events yielding a single neutral H or D. 

After the ion beam is injected into the TSR, phase- 
space cooling and radiative cooling of the molecular ions 
are allowed to proceed over a period of typically 2 s be- 
fore data are taken. The laboratory energy E^^t of the 
electrons in the ETS is then varied in cycles of tens of 
milliseconds between values of E^^t = Ec, which yield 
Ed^Q ("cooling"), oiEe.t = ^r, which yield Ed = 10 eV 
("reference"), and of i?e,t = Em ("measurement"), var- 
ied according to the range of detuning energy E^, to be 
scanned. The rates of DR signals, i?DR, and of single 
mass events, Ri with z = H,D, are recorded at these en- 
ergies. The variation of E„i from injection to injection 
yields the rates as a function of energy, from which to- 
gether with the electron densities the relative DR rate 
coefficient at Em is obtained by 



aDR(-B„ 



^DR(-Em) Ru{Ec) ne^t{Er 



RuiEcEm) RuniEr) ne,t{Em)' 



(3) 



Here, i?DR(£'m) is the DR rate at Em and i?D(-£'c,^m) 
the rate of deuterons at Ec, both averaged over all cycles 
of a given Em, whereas Rd{Ec) and R^niEr) are the 
rates of deuterons and of DR at E^ and E^, respectively, 
averaged over all values of Em of one data set. Differ- 
ent data sets may differ in the ratio of Rd{Ec) / RuRiEr) 
depending on the average pressure in the storage ring. 

As a function of the laboratory energy Em-, the rate 
coefficient aDR measured by this method shows a peak 
at Em = Ec, corresponding to the smallest collision ener- 
gies e where the DR cross section crDR(£) diverges, while 
the rate coefficient aDR(£'m) assumes a finite value de- 
termined by the spread of the electron-ion collision ve- 
locities. From Eq. (pi), relating Ee,t to Ed, it is seen that 
any variations of Ee,t lead to much smaller changes in the 
collision energies, distributed in a narrow range around 
Ed- Hence, sorting the data into suitable intervals of Ed, 
highly resolved variations of the energy resolved merged- 
beams rate coefficient a-oKiEd) are measured. This rate 
coefficient is linked to the DR cross section via 



O^DRiEd) 



va-DR{£)pi£,Ed)de 



(4) 



where v — (2e/?7ie)^^^ is the collision velocity and 
p{e, Ed) the distribution of collision energies for a given 
detuning energy, which follows from the anisotropic ther- 
mal velocity distribution and reflects the electron beam 
temperatures T± and Ty . For the transformation between 
-Ee,t and Ed, small long-term drifts between the voltages 
supplied in the ETS and the electron cooler are moni- 
tored by regular scans over the peak at Em = Ec, fixing 
Ed = at the maximum position of the peak in the 
analysis of the data taken in the temporal vicinity. The 
weighted average of these rate coefficients for the differ- 
ent data sets is then corrected for toroidal effects [33] 
and finally absolutely normalized to the rate coefficient 
aDR(£'d = 9.8 eV) = 6.6 x 10"^ cm^s'^ reported by Al- 
Khalih et al. HI]. 
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FIG. 1. (Color online) DR rate coefficient of HD"^ in the 
range of 10"* to 20 eV, measured with a cold electron beam 
from a photocathode (fcTx =0.6(2) meV, black line and sym- 
bols). This is compared to the previous benchmark measure- 
ment [11] {kT±—2 meV) shown as a violet (gray) line. 



In addition to the energy-resolved velocity-averaged 
rate coefficient aY)^{Ed), also the reduced rate coeffi- 

1 /? 

cient aDR(£^d)£'d is used for the presentation of exper- 
imental and theoretical data. This choice suppresses the 

1/2 

overall decrease of aDR,(£^d), scahng as ex 1/E^' , which 
is caused by the 1/e behaviour of the non-resonant DR 
cross section. For a cross section crDR(£) = K/e with con- 
stant K, this reduced DR rate coefficient for Ed ^ fesTi 
assumes the constant value oi K{2/me)^^^. Near Ed = 0, 

1 /2 

the value of auR{Ed)E^' is lowered by the convolution 
over the electron energy distribution. For asymptotically 
small Ed, it rises as K {2t: Ed /mekBT^Y^^ because of the 
convolution with the electron energy distribution and 
reaches >90% of its final value at Ed > 4:kBT±, which 
corresponds to ^2.5 meV in the present case. Thus, 

1/2 

above a few meV, variations of auR{Ed)E^' only oc- 
cur if the energy dependence of ctdr deviates from 1/e, 
which emphasizes the resonant contributions. 



III. 



EXPERIMENTAL RESULTS 



The measured merged-beams rate coefficient as func- 
tion of the detuning energy is shown in Fig. [T] Tempo- 
ral changes of the energy dependence became insignifi- 
cant for times larger than 2 s after the injection, so that 
we present data from the storage time interval of 2-10 s. 
High-resolution scans were limited to the energy region 
below 0.7 eV, above which the rate coefficient becomes 
very small. The high-energy peaks near 10 and 16 eV 
were scanned with lower resolution and less measure- 
ment time. Here the previously observed structure |11| 
for vibrationally relaxed HD"*" ions is well reproduced. 
Of the previous benchmark measurements presented by 
Al-Khahh et al. [11], that from CRYRING is chosen for 
comparison with the present data as it was obtained with 
the lowest electron temperature so far reported for this 




25 50 75 100 

Detuning energy E^ (meV) 



125 



FIG. 2. (Color online) Reduced DR rate coefficient of HD+ 
in the low energy range, measured with a cold electron beam 
from a photocathode (fcrL=0.6(2) meV, black line and sym- 
bols). This is compared to the previous benchmark measure- 
ment [jA, {kT±=2 meV) shown as violet (gray) line. Newly 
resloved structures are highlighted by arrows. The horizon- 
tal bars indicate the collision energy spread (FWHM) of the 
photocathode measurement in the respective energy regions. 



system, kgT^ = 2 meV. In the present results, the struc- 
ture at low energies appears much better resolved. Thus, 
the peaks at 3, 8, and 19 meV become more pronounced 
and sub-structures appear much more clearly, as illus- 
trated in Fig. [T] by the additional low-energy peak near 
2.5 meV. 

A more detailed view of the low-energy structure is 
presented in Fig. [2] using the reduced rate coefficient 

1 /2 

ctDR{Ed)E^ . This reveals clear sub-structures and 
newly observed peaks between 30 and 100 meV. 

For a determination of the electron temperature, the 
low-energy edges of the lowest- lying peaks, corresponding 
to energy regions near 2, 7, and 17meV, were compared 
to a model rate coefficient. It was obtained from a cross 
section composed of an overall 1/e decrease and three 
infinitely narrow resonances near the peak positions, 
folded with the anisotropic thermal velocity distribution. 
For appropriately chosen peak positions, good agreement 
was found for ksT^ = (0.6 ± 0.2) meV, while the high- 
energy edges are consistent with fcsTJi — 0.025 meV. 
This implies a collision energy spread (full width at 
half maximum, FWHM) of 1 meV or less in the low- 
energy limit, ~2 meV at Ed = lOmeV and ^6.5 me V 
at Ed = 200 meV, scahng as oc e]/^ for Ed > 80meV. 
The energy resolution obtained in the present measure- 
ments is thus comparable to, and for very low collision 
energies even below the energy splitting of the two lowest 
rotational levels of IID+ of 5.44 meV, so that the rota- 
tional structure of the ro-vibrationally excited Rydberg 
resonances ^36^ in the electron collision continuum can be 
expected to become visible. 

After leaving the ion source the HD+ ions are vibra- 
tionally and rotationally excited (rotational and vibra- 
tional initial ionic levels N^ and v^ , respectively). But 
while the vibrations are known to radiatively relax within 
storage times considerably shorter than 1 s [13] , the ro- 



tational excitation of the HD+ ions built on the vibra- 
tional ground state is expected to reach an equilibrium 
only after several seconds of storage, with a number of 
rotational states remaining populated at all times. The 
cooling rate and the resulting population distribution fi- 
nally reached are determined by the interaction with the 
blackbody radiation field of 300 K that pervades the inner 
volume of the room-temperature ion storage ring, and by 
the continuous interaction of the HD+ ions with the elec- 
trons of the electron cooler and the electron target. While 
the Einstein coefficients between the rotational states re- 
quired to describe the interaction of HD+ with the radia- 
tion field are well known [32] , the expected large inelastic 
cross sections for low-energetic electron scattering on ro- 
tational states of HD+ could only recently be verified by 
observing the rotational cooling of IID+ by super-elastic 
electron collisions (SEC) in a dedicated experiment |23) . 
Based on these findings we estimated the time develop- 
ment of the rotational populations by adjusting the calcu- 
lations described by Shafir et al. [53] to the present exper- 
imental situation. In particular, the larger temperature 
of the electron beam of the electron cooler of 16 meV in its 
co-moving reference frame, substantially higher than that 
of the electron beam used in Ref. [23j (2.8 meV), led not 
only to rotational de-excitation by SEC. It also caused 
considerable excitation of the lowest rotational states by 
inelastic collisions with electrons in the high-energetic tail 
of the velocity distribution in the electron cooler. Con- 
versely, cooling and heating of the iV^ distribution by 
inelastic electron collisions in the ETS were found to be 
negligible, mainly because of its much lower electron den- 
sity. While steady-state conditions were found in these 
estimates to be reached only asymptotically at the end 
of the observation window (^10 s), the average N^ pop- 
ulation distribution, weighted by the storage lifetime of 
the HD+ beam of ^7 s, can be approximated by a Boltz- 
mann distribution with a rotational temperature in the 
range of 250-300 K (the range of uncertainty reflecting 
the limited accuracy at which the SEC cross sections 
could be verifled in Ref. [23]). In the following discus- 
sion we therefore assume a thermal distribution of 300 K 
to reasonably represent the relative populations of the 
HD"*" rotational levels during the present measurement. 
The resulting populations of the seven lowest rotational 
states with quantum numbers N^ = 0-6 are compiled 
together with other properties in Table [l] 



TABLE I. Properties of rotational levels in the vibrational 
ground state HD"'", showing A'^^^ (first line), excitation ener- 
gies in meV (second line, from Ref. [34]), relative thermal 
populations at 300 K (third line), and radiative lifetimes in s 
(fourth line, from Ref. |35)). 
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5.44 

0.104 0.251 

140.2 



16.28 
0.271 
14.61 



32.47 
0.199 
4.04 



53.91 
0.108 
1.64 



80.46 
0.046 
0.823 



111.99 

0.015 0.006 

0.469 



IV. 



THEORETICAL METHOD 



Dissociative recombination is a reactive collision in- 
volving ionization channels, describing the scattering of 
an electron on the molecular ion, and dissociation chan- 
nels, accounting for the atom-atom scattering. The main 
steps of our current MQDT treatment [15] are described 
below. 



the nuclear wave-functions corresponding to a dissocia- 
tive state and to an ionization channel, respectively. 

This procedure applies in each A-subspace, and results 
in a block-diagonal global interaction matrix. The block- 
diagonal structure, corresponding to the available A sym- 
metries, propagates to the further built matrices. 



B. Construction of the reaction matrix K. 



A. Construction of the interaction matrix V 

This is performed in the outer shell of the region of 
small electron-ion and nucleus-nucleus distances, that 
is in the 'A-region' |37j . where the Born-Oppenheimer 
context is appropriate for the description of the collision 
system. The good quantum numbers in this region arc 
N, M , and A, associated respectively to the total angu- 
lar momentum and its projections on the z-axes of the 
laboratory-fixed and of the molecule- fixed frames. 

Within a quasi- diabatic representation |361 1381 [55] . the 
relevant states are organized in channels, according to the 
type of fragmentation which they are meant to describe. 
An ionization channel is built starting from the ground 
electronic state of the ion and one of its ro- vibrational 
levels N^v^, and is completed by gathering all the mono- 
electronic states of a given orbital quantum number /. 
These mono-electronic states describe, with respect to the 
iV+u+ threshold, either a "free^' electron - in which case 
the total state corresponds to (auto) ionization - or to a 
bound electron - in which case the total state corresponds 
to a temporary capture into a Rydberg state. In the A- 
region, these states may be modeled reasonably well with 
respect to the hydrogenic states in terms of the quantum 
defect ii{^{R), dependent on the internuclear distance R, 
but assumed to be independent of energy. 

An ionization channel is coupled to a dissociation one, 
labeled dj, by the electrostatic interaction l/ri2. In 
the molecular-orbital (MO) picture, the states corre- 
sponding to the coupled channels must differ by at least 
two orbitals, the dissociative states being doubly-excited 
for the present case. We account for this coupling at 
the electronic level first, through an i?-dependent scaled 
"Rydberg-valence" interaction term, Vj'^j , assumed to 
be independent of the energy of the electronic states per- 
taining to the ionization channel. Subsequently, the inte- 
gration of this electronic interaction on the internuclear 
motion results in the elements of the interaction matrix 
V, 

V^^%,,AE,E) = (xL,l^i;,l^lxf.„.), (5) 

where E is the total energy and Xjvd ^^'^ Xn+v+ ^^^ 



Starting from the interaction matrix V, we build 
the reaction K-matrix, which satisfies the Lippmann- 
Schwinger integral equation 



/c = v + v 



1 



E-Hn 



K. 



(6) 



This equation has to be solved once V - whose elements 
are given by Eq. ([5]) - is determined. Here Hq is the zero 
order Hamiltonian associated to the molecular system, 
while neglecting the interaction potential V. It has been 
proven that, provided the electronic couplings are energy- 
independent, a perturbative solution of Eq. ^ is exact 
to second order 140 . 



C. Diagonalization of the reaction matrix K, 

In order to express the result of the short-range in- 
teraction in terms of phase-shifts, we perform a unitary 
transformation of our initial basis into eigenstates. The 
columns of the corresponding transformation matrix U 
are the eigenvectors of the K-matrix, and its eigenvalues 
serve to evaluate the non-vanishing elements of a diago- 
nal matrix rj: 



K\J 



tan(T7)U. 



(7) 



D. Frame transformation to the external region 

In the external zone - the 'B-region' [37', characterized 
by large electron-core distances - the Born-Oppenheimer 
model is no longer valid for the whole molecule, but only 
for the ionic core. A is no longer a good quantum number, 
and a frame transformation |41H43| is performed between 
coupling schemes corresponding to the incident electron 
being decoupled from the core electrons (external region) 
or coupled to them (internal region). The frame trans- 
formation coefficients involve angular coupling, electronic 
and vibrational factors, and are given by 



■'lN+v+,Aa 



2N+ + 1 



1/2 



2iV+l 

l + r+T(-l) 






{1{A-A+)N+A+\IN+NA) 



(8) 



(9) 



as well as by Si]y+y+\a ^nd Sd^, Aai which are obtained 
by replacing cosine with sine in Eqs. (|8| and Q. In 
the preceding formulae, Xn+v+ ^^ ^ vibrational wavefunc- 
tion of the molecular ion, and Xnv i^ ^ vibrational wave- 
function of the neutral system adapted to the interaction 
(A) region. The matrix elements (x^+^+| cos(7r^^(i?) + 
t]^)\Xnv) represent the ro-vibronic couplings giving rise 
to the resonant capture processes. The quantities t'^ 
and T are related to the reflection symmetry of the ion 
and neutral wave function respectively, and take the val- 
ues +1 for symmetric states and —1 for antisymmetric 
ones. The ratio in front of the sum in the right mem- 
ber of Eq. (l8| contains the selection rules for the rota- 
tional quantum numbers. The indices dj (j = 0, 1, 2, . . .) 
stand for the states of a given symmetry, open to dissoci- 
ation at the current energy, a denotes the eigenchannels 
built through the diagonalization of the reaction matrix 
K and — tan(r7^)/7r, U^ ^ are its eigenvalues and the 
components of its eigenvectors respectively. 

The projection coefficients shown in Eqs. (|8| and ^ 
include the two types of couplings controlling the pro- 
cess: the electronic coupling, expressed by the elements 
of the matrices U and t], and the non-adiabatic coupling 
between the ionization channels, expressed by the matrix 
elements involving the quantum defect iif'{R). This lat- 
ter interaction is favored by the variation of the quantum 
defect with the internuclear distance R. 



E. Construction of the generalized scattering 
matrix X 



The matrices C and S with the elements given by Eqs. 
((sl) and ^ are the building blocks of the 'generalized' 
scattering matrix X: 



X = 



C + iS 
C-iS' 



It involves all the channels, open ("o") and closed ("c"). 
Although, technically speaking, the open channels only 
are relevant for a complete collision event, the partici- 
pation of the closed channels may influence strongly the 
cross section, as shown below. 

The X matrix relies on 4 block sub-matrices, 



X = 



Xqo Xq 
Xrn X^ 



(11) 



I 

where "o" and "c" label the lines and columns corre- 
sponding to open and closed channels, respectively. 



F. Elimination of closed channels 

The building of the X matrix is performed indepen- 
dently on any account of the asymptotic behaviour of 
the different channel wavefunctions. Eventually, impos- 
ing physical boundary conditions leads to the 'physical' 
scattering matrix, restricted to the open channels f44j : 



1 



"^cc — cxp(— i27rf) 



(12) 



It is obtained from the sub-matrices of X appearing in 
Eq. ( |11[ ) and on a further diagonal matrix v formed with 
the effective quantum numbers i^jv+i>+ = [2(£'Ar+u+ — 
_E)]^^/^ (in atomic units) associated with each vibra- 
tional threshold Ejq+y+ of the ion situated above the 
current energy E (and consequently labelling a closed 
channel) . 



G. Cross section evaluation 



For a molecular ion initially in the level N^v^ and 
recombining with an electron of kinetic (collision) energy 
e, the cross section of capture into all the dissociative 
states dj of the same symmetry is given by 



A^, sym 



2iV+l 



sym 



4£2iy^ 



1 






NA 



(13) 



On the other hand, the cross section for a ro- vibrational 



transition to the final level Ntvf, giving coUisional 



(10) (de-)excitation, is 



N, sym 


TT 2N + 1 

" 4e 2N+ + 1 







syin 



(14) 



Here p^^™ is the ratio between the multiplicities of the 
neutral and the target ion. After performing the MQDT 
calculation for all accessible total rotational quantum 
numbers N and for all relevant symmetries, one has to 



add up the corresponding cross sections in order to ob- 
tain the global cross section for dissociative recombina- 
tion or ro- vibrational (de-)excitation as a function of the 
electron collision energy e. 



V. COMPUTATION OF CROSS SECTIONS 
AND RESONANCE ASSIGNMENT 

A. Computations 

Since the lowest ^I]+ doubly excited state has by far 
the most favorable crossing with the ion curve for colli- 
sions taking place at low energy, our analysis has focused 
on its contribution to the DR cross section. To improve 
the accuracy of the cross sections, we have used for this 
symmetry the adjusted quasi-diabatic molecular data - 
electronic potential curves and couplings - as described 
in Ref. |15j . In particular, the quantum defect functions 
were obtained by the procedure of Ross and Jungen [IS] 
using the large set of adiabatic potential curves from Det- 
mer et al. |46j . With these data, we have performed 
MQDT-based [IHl IMl US] computations of the DR cross 
section for vibrationally relaxed HD+ on its first 13 ro- 
tational levels {N^ = 0-12, ■y^ = 0), a range which is by 
far sufficient to satisfactorily model an equilibrium rota- 
tional distribution corresponding to temperatures below 
1000 K. 

For each initial rotational state N^ of the ground 



„+ 



— 0, the cross section a 



diss^JV; 



vibrational level 

has been obtained by summing up the contributions 

a '^^"^ J. over all accessible values of N and all sym- 

metrics. Once the total cross-sections were computed, 
we have calculated Boltzmann averages for different ro- 
tational temperatures between 100 and 1000 K. This out- 
put has then been used for a further average, consisting 
of the convolution with the anisotropic Maxwell distri- 
bution of the relative velocities between electrons and 
cations in storage ring conditions. This yields the energy- 
resolved rate coefficient aDR(-Ed) described by Eq. Q 
suitable for comparison with the experiment, where ap- 
propriate electron temperatures T± and T^\ must be set 
to define the collision energy distribution p{e,Ej). For 
the calculations presented here we have chosen /cbTl = 
0.5 meV and fesTn = 0.02 meV. Regarding the rotational 
temperature assumed in the Boltzmann average to deter- 
mine the rotational level populations, we will assume the 
value of 300 K considered as a reasonable estimate of the 



experimental conditions (see Sec. Ill I 



The high experimental resolution and the finite num- 
ber of rotational levels giving significant contributions 
to the rate coefficient aDR(i?d) suggest to investigate the 
origin of individual structures in the theoretical cross sec- 
tion and their correspondence to the energy dependence 
of the experimental cross section. Such structure can 
in particular be of resonant character, representing elec- 
tronically bound Rydberg states in the various ionization 



channels as defined above, which are associated to well- 
defined quantum numbers iV+ , w"*" of the ion core in the 
'A-region'. Single rotational states N^ of the incident 
ion are discussed separately before the average over the 
rotational level populations is added as a final step. The 
choice of the electronic orbital momentum / of the Ryd- 
berg resonances corresponds to that applied earlier |20) . 
where it is restricted to / = and 2 (s and d states, only). 
The resulting angular momenta N that contribute to the 
sum in Eq. ( 13 1 as well as the channel quantum numbers 



TV"*" are those given in Table 1 of Ref. 



B. Resonance assignment 

The assignment procedure has been carried out for 
each of the first six rotational states N^ of the ground 
vibrational state of HD+(2i;+), that is, 7V+ = 0-5. As 
suggested by Eq. Q, the energy dependent DR cross 
section o'DR(e) — '^^■'^^™,^+ +{^) is represented in the 

form of W(TDR(e), which is denoted as the "local" rate 
coefficient. The step to the experimental energy-resolved 
rate coefficient then consists only in the convolution with 
p{e,Ed) to obtain aY:iK{Ed)- For simplicity, the same 
energy scale (denoted as "electron collision energy") is 
used in the figures both for the collision energy e itself 
and for the detuning energy Ed scanned in the experi- 
ment as a proxy for e; moreover, the convolved energy- 
resolved rate coefficient aYy^{Ed) is represented on the 
same vertical scale as the unconvolved local rate coeffi- 
cient i'crDR(e). The fully assigned local rate coefficients 
and the results of the convolution are presented in Fig. [3] 
and 4 for HD"'"(^S+,v^ = 0) in the initial ro-vibrational 
states N^ = 1 and 2, respectively, which have the highest 
relative populations at 300 K (see Table ll| . The energy 
range of 0-24 meV is chosen in order to focus on the reso- 
nances that determine the DR rate coefficient under cold 
plasma conditions up to ^100 K. 



1. Results for ttie most populated initial rotational states 

In Fig. [3] the calculated energy-resolved DR rate coef- 
ficient is presented for NJ^ — 1. The result aDR(^d) of 
the convolution in Fig. ^ a) shows a number of maxima 
and minima whose assignment has been obtained by con- 
sidering the unconvolved local rate coefficient vaunis) in 
Fig.lSlb) and, at the most elementary level shown in Fig. 
IsFc), the contributions to it arising from the accessible 
states of total angular momentum N. Details of the res- 
onance assignment will be discussed in Sec. |VB 2| It is 
guided by considering the Rydberg series corresponding 
to various relevant ionic-core excitations as indicated by 
the prediction bars in the lower part of the figure. Purely 
rotational excitations from N^ = 1 to iV+ = 3,5 lead to 
high-n Rydberg resonances in the near-threshold elec- 
tronic continuum. Only small contributions to auniEd) 
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FIG. 3. (Color online) Calculated and assigned MQDT rate 
coefficient for the DR of HD"*" (^ E^ , ii^ — 0) in the initial state 
A^,^ = 1. (a) Convoluted rate coefficient a-DB.{Ed) as calcu- 
lated (upper thick curve) and multiplied by the expected rel- 
ative Boltzmann population in A''^^ = 1 (lower thick curve) 
compared to the measured data (dots and thin line). Hori- 
zontal bars: FWHM experimental collision energy spread (see 
Sec. nil. (b) Local rate coefficient vaoR{s) summed over the 



contributions of all total angular momenta A'^ (full line) and its 
convolution (dashed line); the same ucrDR(e) is also shown in 
(a) by a blue (light gray) curve, (c) Contributions to vaonie) 
from TV = 1-3. At the bottom the energies e corresponding 
to Rydberg states n of the relevant channels with core quan- 



tum numbers A'^ 



are indicated by black vertical bars for 



d states and red (gray) bars for s states (numbers for n) with 
the energies from Eq. ( 15 1 and ( 16 1. Assigned resonance min- 
ima and maxima are labeled hyjN"^ , v^ , n) . Normal type is 
used for d and italic type for s Rydberg states. 



result from N^ = 3 (n > 23), but significant effects ensue 
for iV+ = 5, n = 14, 15 even after the convolution. Signif- 
icant peaks and dips in the rate coefficient are also caused 
by higher ro-vibrational core excitations, which lead to 
Rydberg resonances with lower principal quantum num- 
ber n. These structures are considerably shifted in energy 
from their positions predicted by the Rydberg formula 
and need more detailed verifications in order to identify 
their origin (see Sec. VB2). Important structures can 
be assigned to the resonances 8d and 4d in the excitation 



channels 



,+ - 



1,A^+ = 1 and v+ = 4,N+ = 1, 



respec- 



tively, whose energies are well known [9] to lie in the 
electronic continuum relevant for the low-energy DR of 
HD+. Note that these cases do not involve any change in 
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FIG. 4. (Color online) Calculated and assigned MQDT rate 
coefficient for the DR of HD^(^E^, vf — 0) in the initial state 
N^ — 2. See caption of Fig. Islfor the notation. 



the rotation, as A^+ = Nlj^ . Moreover, the 3d resonance 
in the even higher «+ = 7, N'^ = 3 excitation channel, 
energetically predicted below the ionization threshold for 
this channel, is found to cause a strong dip in the rate co- 
efficient at only ~2 meV collision energy for this incident 
rotational level. 

The energy-resolved DR rate coefficient for A^^^ = 2, 
and the contributions to it of the individual total angu- 
lar momenta N, are shown in Fig. |4] Considering the 
core excitations indicated in the lower part of the plot, 
again some resonances from purely rotational excitation 
remain visible even after the convolution over the elec- 
tron energy distribution. Thus, the 13d Rydberg level of 
the excitation channel Nl^ = 2 — )• N^ = 6 leads to a 
clear structure around 16meV; in addition, the 19s and 
19d Rydberg levels of the channel N^ = 2 -)■ N+ = A\ie 
very close to the 1{D'^{N^ — 2) ionization threshold and 
likely cause the calculated aBR{Ed) to rise by a factor of 
>4 in the range of Ed < 0.5 meV. The dominant struc- 
tures in the predicted aj:)B.iEd) are again the 8d and Ad 
resonances in the channels v'^ = I and v^ ~ 4, respec- 
tively, with N+ = N:;^ = 2. The n = 3 states of the 
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FIG. 5. (Color online) Calculated A'^ = 1 contribution 
to the MQDT local rate coefficient u(tdr(£) for the DR of 
HD+(^E+,u+ = 0) in iV+ = 1 and illustration of the valida- 
tion for the assigment of low-n resonances. Black line: full 
result for A^ = 1; red (gray) line: result obtained on eliminat- 
ing the excitation channel {N'^,v'^) — (1,1) as symbolically 
indicated in the right margin. The resonances are labeled 
using the notation of Fig. [S] 



channels w+ = 7, A'^^ = 4 and v~^ = 8, N'^ = predicted 
in this energy region probably lead to an additional broad 
modulation of the calculated local rate coefficient. 

The comparisons to the experiment in Fig. pTa) and 
Fig. Ka) show that the calculations for only the two most 
populated initial rotational states already well reproduce 
the observed strong structure near 20meV, which coin- 
cides with the predicted Sd, w+ — 1 peaks for both Nl^ . 
The calculated positions of these peaks are fairly inde- 
pendent of N^ as the related resonant capture processes 
do not involve a change in the core rotation and n is suf- 
ficiently high, so that variations of the quantum defect 
with N'^ have only a small effect. This is different for 
most other resonances. The calculated rate coefficient 
after summation over N^ will be presented in Sec. 
where also the main resonant contributions for the otT 
initial rotational states will be labeled. 



2. Procedure and verification of resonance assigment 

The resonance assignment demonstrated in the pre- 
vious section is the result of a two-step strategy, based 
on the prediction of the resonance energies and a valida- 
tion of the resonant effect onto the local rate coefficient. 
These two steps are implemented as follows. 

(i) Prediction. Considering a channel with the core 
quantum numbers N~^ and w"*", the conservation of the 
total energy in the resonant electron capture process into 
a Rydberg state with principal quantum number n and 
orbital angular momentum I can be approximated by 



e + E{N+,v: 



E N^ 



n 



{n- hn+v+iY 



(15) 



Here TZ is the Rydberg energy constant, flpf+y+i the suit- 
ably averaged quantum defect of the relevant channel (see 
below), and e the energy of the incident electron associ- 
ated with the resonance nl of channel TV"*- , v~^ . These 
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FIG. 6. (Color online) Same rate coefficient as in Fig. [5] 
(black line) and validation tests, shown as red (gray) lines, by 
(a) eliminating the (1,4) excitation channel and (b) slightly 
enlarging the quantum defect of the d Rydberg levels. Reso- 
nance labels use the notation of Fig. IS] 



energies are used to create the Rydberg series prediction 
bars in the figures of this Section. 

The quantum defect can in lowest approximation, ap- 
propriate for low v~^ and high n, be calculated at the 
equilibrium internuclear distance R of the ground-state 
molecular ion. However, for lower n states, which be- 
come relevant at higher w+, the anharmonicity of the 
HD"*- potential curve together with the i?-dependence of 
the quantum defects leads to imprecise predictions of the 
Rydberg resonance energies. Hence, we obtain the quan- 
tum defects for a particular excited core state 7V+, -y"*- as 
the average 



IJ-N+v+l — iXN+v+lfJ-l {R)\XN+v+)- 



(16) 



The actual structure of the local rate coefficient due to 
the Rydberg resonance is still shifted from the predic- 
tion of Eq. ( 15 ) through the perturbation induced by the 
dissociative state, which is properly taken into account 
by the MQDT formalism. Hence, in a number of cases 
structures cannot be uniquely identified using the predic- 
tion bars alone, and the assignment to specific excitation 
channels and Rydberg resonances requires a validation 
procedure. 

(ii) Validation. For a unique identification of the ex- 
citation threshold N^,v^, we eliminate the respective 
resonant contributions in the cross-section calculation. 
For this purpose, the matrix elements representing the 
ro-vibronic couplings for this specific threshold are set to 
zero in Eq. ^ for both C;Ar+«+,Aa and 5/jv+i,+,Aq- In 
most cases the comparison between the cross section ob- 
tained in these trials and the original calculated cross 
section enables us to identify the structures due to the 
considered threshold. Similarly, the contributions from 
different resonance orbital angular momenta I are distin- 
guished by temporarily setting the corresponding quan- 
tum defect function fJ.f{R) to zero in Eq. (ISl. 

A further tool for resonance identification is to probe 
the sensitivity of the energy at which a structure occurs 
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FIG. 7. (Color online) Same rate coefficient as in Fig.ls] but 
for A'^ = 3 (black line) and performing validation tests, shown 
as red (gray) lines, by (a) eliminating the (3, 7) excitation 
channel, (b) eliminating the resonant capture into s Rydberg 
levels, and (c) reducing the quantum defect of the d Rydberg 
levels. Resonance labels use the notation of Fig. [3] 



to slight changes Aflj^+y+i of the quantum defect, which 
lead to corresponding changes of the resonance energy 
according to 



Ae 



27^ 



{n 



fJ-N+v + l)^ 



A/2 



N+v + l- 



(17) 



Low-lying mono-excited Rydberg states are strongly cou- 
pled with the dissociation continuum, which results in 
large width of the resonance profile. Furthermore, the 
inclusion of the vibronic coupling via the frame trans- 
formation induces asymmetric resonance profiles. Thus, 
the assignment of such states using only predicted reso- 
nance positions is normally not possible, and the detailed 
validation steps discussed here are needed. 

The validation of the peak and valley positions for 
the 8d resonance in the channel {N^,v^) — (1,1) for 
N^ = 1 is demonstrated in Fig. pi The threshold elimi- 
nation procedure is used on the (1, 1) excitation channel 
and leads to a disappearence of the presumed {l,l)8d 
minima and maxima. The corresponding validation of 
the (1, 4)4d minimum for the same initial state is demon- 
strated in Fig. l6ja) regarding the associated excitation 
channel, while Fig. WCh) shows the effect of a shift of the 
quantum defect for the d Rydberg states, which leads to 
downward displacements of both the presumed (l,4)4d 
and (1, l)8d structures, with a larger shift observed as 
expected for the structure assigned to Ad. 

All three steps were combined to identify the reso- 
nance (3, 7)3d in the iV = 3 contribution to the N^ — 1 
rate coefficient, as shown in Fig. [7] Near 1.2 meV, the 
MQDT local rate coefficient shows more mimima than 
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FIG. 8. (Color online) Calculated N — 2 contribution to 
the MQDT local rate coefficient for DR of HD+(^E+,i'+ = 
0, N^ = 0) ions, (a) Result with resonance labels using the 
notation of Fig. |3] The lower traces show the validation pro- 
cedure of the resonance assignment in particular for the mini- 
mum and the maximum marked by the rectangle, which turn 
out to arise from Rydberg states of different series. Modified 
rate coefficients are shown by red (gray) lines for elimiation 
of the excitation channels {N^,v'^) = (0,8) (b) and (2,4) (c) 
as well as for a reduction of the d quantum defect (d). The 
displaced rectangle in (d) marks the shifted (0, 8)3d minimum. 



energetically predicted at these energies (only two min- 
ima should occur along the (3, 0) Rydberg series from 
23d and 23s). The third minimum could arise from the 
nearby (3, 7)3s resonance or possibly the (3, 7)3d reso- 
nance occuring at e < 0. On removal of all s reso- 
nances in Fig. I7[b), only the narrow structure close to 
the predicted (3^0) 23s energy disappears. The deepest 
minimum near 1.2 meV, occurring between (3, 0)23d and 
(3, 0)23s, clearly has d character as it remains when the s 
contribution is removed. Also, on slightly lowering the d 
quantum defect in Fig.m c), this minimum moves upward 
by a larger amount than the nearby (l,4)4d, identifying 
it as (3, 7)3d. 

As a final example for the validation procedure, the 
identification of resonance maxima and minima is demon- 
strated in Fig. [8] for the MQDT local rate coefficient 
for the DR of HD+ ions in the rotational ground state 
{N^ = 0). In addition to the prediction bars shown, 
the (0, 8)3d and (2,4)4(i states are calculated from Eqs. 
(15) and ( [l6| to he out of scale, at e = — 4.49meV and 
e = 24.96 mcV, respectively. Elimination of the (0, 8) 
excitation channel in Fig. [slb) leads to disappearence 



of the minimum near 19meV, while the elimination of 
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FIG. 9. (Color online) Reduced rate coefficient for the DR 
of RD+{^T,j,v+ = 0). Theory after the convolution of Eq. 
with kT± = 0.5 meV and fcTy = 0.02 meV. (a) Contri- 
butions of the six lowest rotational levels, weighted by the 
thermal Boltzmann populations at 300 K (broken lines with 
numerical labels for N^), total MQDT rate coefficient (thick 
full curve), and experimental DR rate coefficient (grey dots 
and thin line). Correspondence between experimental and 
theoretical peaks is indicated by gray connecting lines, (b) 
Assignment of the dominant contributions to the peaks in 
the A'^j^-summed MQDT rate coefficient, using the notation 



NtiN' 



with normal type for d and italic type for s 



Rydberg states. Blue (light gray) curve: calculated total rate 
coefficient «(TDR(e)e^ before convolution. 



the (2, 4) excitation channel in Fig. Isjc) suppresses the 
nearby maximum. On a shght decrease of the d quantum 
defect in Fig.^d), both the minimum and the maximum 
are upshifted, revealing their d character; as expected be- 
cause of the lower principal quantum number involved, 
the minimum assigned to (0, 8)3d shifts more strongly 
than the maximum assigned to (2, 4)4d. 



VI. DETAILED COMPARISON BETWEEN 
STORAGE-RING AND MQDT RESULTS 

For the comparison between experimental and theo- 
retical results, we choose the reduced DR rate coefficient 

1/2 

ctDR{Ed)E^ which emphasizes the resonant variations. 



The Boltzmann-weighted contributions of the individual 
HD+ rotational states as well as the total MQDT rate 
coefficient obtained by summing these contributions are 
shown in Fig. ^a), while the main resonances assigned 
to the structures in the total MQDT rate coefficient are 
indicated in Fig. Wih). For the shape of the theoretical 
rate coefficient, the N^ = 1 contribution turns out to 
have the largest influence. In fact, the predicted dou- 
ble structure between 4 and 7meV resulting from the 
(l,4)4(i and (5, 0)14(i resonances for this N^ remains 
visible also in the sum over all initial rotational states. 
HD"*" ions in N^ = 1 are also predicted to yield the 
largest contribution to the peak near 19meV through 
the (1, l)8d resonance. This resonance is further empha- 
sized by other resonances of the same Rydberg level and 
with 7V+ = N^ , which occur at essentially the same e for 
N^ = 1-4. Near 14 meV the maxima and minima caused 
by the Rydberg levels (5, 0)15d of iV:+ = 1 and (6, 0)13d 
of N^ = 2 add up. Theory also predicts two structures 
at very low energy s < 1 mcV: a peak due to combined 
ro-vibrational excitation of N^ = 3 ions, caused by a 
7d Rydberg level, and a rather strong 16d Rydberg res- 
onance through purely rotational excitation of N^ — 
ions to N^ = 4. The latter prediction would significantly 
influence the DR of HD+ in a low-temperature environ- 
ment, where iV^^ = will be much more important than 
in the present case. [Note that also the n = 29 Rydberg 
levels of the N^ — 2 excitation channel occur in very 
close vicinity of the ionization threshold (e = 0) for these 
ions, as shown in Fig. [8j but the MQDT calculation pre- 
dicts lower contributions to the DR rate for these states.] 
Altogether, while the strongest resonances in the calcu- 
lated rate coefficient are due to combined ro-vibrational 
excitation, also purely rotational excitation and capture 
in Rydberg levels with n = 13-16 are found to give sig- 
nificant contributions. 

In comparison to the experiment, a clear correspon- 
dence is observed between the calculated peaks of N^ — 
1-4 near 19meV and the dominant experimental struc- 
ture found at this energy. At lower energies, less good 
agreement is found; in particular, the prominent experi- 
mental peaks at 2-4 meV (with a sub-structure) and at 
8 meV are not correctly reproduced by the theory. Con- 
sidering the trend of low-n resonances to acquire a large 
width, it appears not unreasonable to suggest a corre- 
spondence between the strong theoretical (l,4)4(i reso- 
nance for N^ = 1 at 6 meV with the experimental struc- 
ture observed 2 meV higher in energy. It can also be ra- 
tionalized that for the low n involved in such resonances, 
both energetic position and width will sensitively depend 
on the underlying quantum defect of this channel and on 
its residual n dependence. On the other hand, this as- 
sumes that the lower experimental peak (or unresolved 
peaks) are caused by levels not properly accounted for 
in the calculation. It must be noted that for a predicted 
high-n resonance, such as (5,0)14(i for iV^^ = 1 near 5 
meV, variations of the quantum defect will not signifi- 
cantly change the resonance position on the considered 
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FIG. 10. (Color online) Reduced DR rate coefficients for the DR of HD"''(^E^,i;^ = 0) from experiment (black dots and thin 
line) and MQDT theory (thick solid curve) after collision energy convolution and initial state averaging as in Fig. [9] Prediction 
bars: Rydberg resonance energies (labeled by principal quantum numbers n; zero quantum defect) below the energy thresholds 
for vibrational excitation to levels v'^ as indicated, assuming AA''^ = in the resonance formation. Near the «+ = 1 threshold, 
its energetic shifts due to molecular rotation are indicated for initial states A^'^^ = 2, 3 and the cases AA'^"'" = and AA^^ — —2. 



Horizontal bars: FWHM experimental collision energy spread in the respective energy regions (see Sec. Ill I 



energy scale (even though the resonance strength and 
the exact position of the associated minima and max- 
ima will be affected). Hence, the 2-4 meV experimental 
peak can be plausibly assumed to be partly caused by the 
(5, 0)14(i, N^ — 1 resonance, while other contributions 
to it are obviously missed by the calculation. 

The theoretical results obtained in the same approach 
over a broader energy range, including the excitation 
limit to the first HD+ vibrational level at 0.237 eV, are 



to vibrational excitation connected with rotational de- 

j+ = 1 threshold 



shown in Fig. 10 and compared to the experimental data, 
using again the reduced rate coefficient aY)R{Ed)E^^ . 
It should be noted that the thermal rate coefficient of 

- Q .S — 1 



ath 



(9 ± 2) X 10^'-* cm^s-i X (300K/r)i/2 deter- 



mined in the limit of low electron temperature by Al- 
Khalili et al. [llj translates into a constant low-energy 

limit of aDR{Ed)Ey^ = {irkBTy^'^ath/^ that amounts 
to (0.13 ± 0.03) X IQ-^^cm^s-^eVi/^, in good agree- 
ment with the average value of the reduced rate coef- 
ficient observed at <0.02eV, while the overall trend of 

1/2 

aDR{Ed)E/ then shows a lower average up to ^0.2 eV. 
A further decrease occurs as expected at the first vibra- 
tional excitation threshold, when the autoionization into 
the electronic continuum of HD+(w+ — 1) becomes ener- 
getically allowed following resonant electron capture (res- 
onant vibrational excitation). 

The energetic thresholds for autoionization plotted in 
Fig. [TOJ which are complemented by the attached Ryd- 
berg series of capture resonances, assume A^^ = inci- 
dent ions and are slightly shifted by rotational effects. 
Thus, for N^ > and autoionization into the continuum 
with A^+ = N^ , the reduction of the vibrational con- 
stant by centrifugal effects moves the threshold positions 
downward by a few meV, resulting in shifts of —0.74 and 
-1.48 me V for Af + = 2 and Af + = 3 incident HD+ ions, 
respectively. The shifts become larger for autoionization 
into the continuum with A^^ = A^^ — 2, corresponding 



where they amount 
Since processes 



excitation and are illustrated for the 
and the case AA^+ = —2 in Fig. 
to —16.3 and —27.3 meV, respectively, 
with AA^+ y^ are clearly important in the MQDT cal- 
culations, as discussed in detail in Sec. |V] all resonance 
positions marked here (in contrast to Sec. Wj can have 
only indicative character considering the rotational exci- 
tation relevant in the present study. With this in mind, 
for Fig. [To] also the quantum defects were neglected in the 
Rydberg state positions marked which, however, mainly 
affects the low-n resonances for w+ = 4 and 8. 



The present MQDT result yields overall a remarkable 
agreement with both the magnitude and the structure 
of the observed rate coefficient. In fact, up to ~0.2eV 
most prominent structures of the reduced rate coefficient 
are resonably well reproduced by the theory with lo- 
cal energy shifts below 5meV. Also an intermittant dip 

1 /? 
in ar,B,{Eci)E^' between 0.10 and 0.14 eV is reasonably 

borne out by the MQDT prediction. Both in experiment 
and theory, the drop in the reduced rate coefficient at 
>0.2eV, i.e., already below the rotationless vibrational 
excitation threshold, is remarkable. In fact, the onset 
of this decrease occurs just at the expected threshold 
energies of AA^+ = — 2 processes for A^^^ = 2 and 3 in- 
dicated in Fig. [TOJ which suggests that resonant electron 
collisions followed by autoionization in this energy region 
lead to vibrational excitation together with rotational de- 
excitation. The initial rotational levels for which this ef- 
fect can be expected are among the most populated ones 
at 300 K (see Table m. 
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VII. DISCUSSION AND OUTLOOK 

At the high energy resolution and theoretical accu- 
racy achieved in the present work, dissociative recombi- 
nation experiments using a stored beam of vibrationally 
relaxed HD+ ions turn out to be a sensitive tool for 
studying rotational effects in resonant low-energy colli- 
sions with molecular cations. The energy-resolved rate 
coefficient reveals a rich structure from ro-vibrationally 
excited Rydberg resonances which drive the DR processes 
observed here as well as the other collision channels of 
ro-vibrational (de-)excitation and elastic scattering. Ex- 
tensive MQDT calculations yield a high level of agree- 
ment with the observed structures, which confirms the 
good quality of the molecular data applied, the efficiency 
of the theoretical approach for predicting the low-energy 
resonant processes, and the assumptions about the rota- 
tional level population underlying the interpretation of 
the experimental results. 

The rotational level population, close to the thermal 
equilibrium caused by radiative transitions in the black- 
body radiation field of the 300-K (room-temperature) 
storage-ring and by the inelastic and superelastic elec- 
tron collisions in the electron cooler, leads to ~5 similarly 
occupied HD+ initial levels with excitation energies well 
above the experimental energy resolution. The largest 
structure, a narrow resonance near 19 meV, is assigned to 
the rotationally excited initial states N^ — 1-4, whereas 
no contribution from N^ = is predicted for this peak. 
The molecular rotation N^ is conserved in the formation 
of these electron capture resonances. The calculation pre- 
dicts similarly large rate coefficients for processes with 
AiV+ — ±2, ±4 which, however, for the various initial 
rotational levels are strongly dispersed over the energy 
scale studied here. While the experimental structures 
are not reproduced in all their detail, the agreement be- 
tween the measured structures and those in the MQDT 
rate coefficient after convolution with the experimental 
energy distribution is remarkable and by far the best for 
any molecular ion to-date, especially at this resolution. 
Not only can most observed structures between 5 and 
200 meV be individually assigned to calculated ones, but 
also the ratio of the calculated and the measured rate 
coefficients remains within 0.5 to 2.0 over essentially this 
entire energy range. 

The ro-vibrational resonant structure in the low-energy 
DR rate coefficient of HD"*" is among the most pro- 
nounced ones observed for this process on any molecu- 
lar cation. Since ro-vibrational resonances of the type 
considered here occur for a wide range of diatomic and 
polyatomic ions, their statistically averaged DR rate co- 
efficients in cold media should be expected to be signifi- 
cantly influenced by the magnitudes and the temperature 
dependences of this resonant DR process. However, most 
experiments so far observed rather structureless energy 
dependences in low-energy DR studies of heavier systems 
[T]. It is likely that this can be attributed to the much 
smaller energetic differences between the rotational lev- 



els for these species as compared to the hydrogen cation, 
which leads to an even larger number of rotational lev- 
els populated at the rotational temperatures of at least 
300 K in the collision-energy resolved DR measurements 
so far performed on such systems. This stands in con- 
trast to the low-temperature conditions of many media in 
which DR is important, and much more structure in the 
energy dependence of the DR rate coefficient may well be 
revealed also in these cases with improved control over 
the rotational populations. For HD+, the present re- 
sults predict that the collision energy dependence of the 
DR rate coefficient for rotationally cold {N^ = 0) ions 
will be dramatically different from that observed here 
and in previous experiments with similar, near room- 
temperature, rotational populations. The development 
of a cryogenic electrostatic ion storage ring suitable for 
electron-ion collision experiments with merged beams is 
in progress [47l [48] and will make energy-resolved DR 
studies possible under conditions where rotational exci- 
tation by the thermal blackbody radiation is largely elim- 
inated and the parent ions are rotationally cold. 

On the theoretical side, the present work demon- 
strates the power of the MQDT method for predicting 
low-energy collision cross sections between electrons and 
molecular cations determined by numerous Rydberg cap- 
ture resonances, including a detailed account of rota- 
tional effects. The magnitude and energetic dependence 
of the DR rate coefficient can be well understood on 
the basis of these calculations after averaging over the 
electron energy distribution and the initial HD+ rota- 
tional distribution. Methods for identifying the individ- 
ual resonant contributions in the theoretical cross section 
even in cases of overlapping and strongly broadened res- 
onances have been implemented. This will be of great 
help in future efforts to further improve the quantum 
defect functions and their implementation in the reso- 
nant cross-section calculation. Full experimental control 
over the initial rotational excitation of the HD+ ions, 
in connection with MQDT calculations using the proce- 
dures presented here, can be expected to yield a com- 
plete understanding of the state-to-state resonant colli- 
sion mechanism for this benchmark reaction also in the 
low-temperature limit which is not yet finally explored. 
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